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Abstract. We present the first analysis of an EUV selected sample of hot DA white dwarfs using a new type of 
atmospheric models. These models take into account the interplay between gravitational settling and radiative ac- 
celeration to predict the chemical stratification from an equilibrium between the two forces while self-consistently 
solving for the atmospheric structure. In contrast to atmospheric models with the assumption of chemical ho- 
mogeneity, the number of free parameters in the new models is reduced to the effective temperature and surface 
gravity alone. The overall good reproduction of observed EUV spectra reveals that these models are able to de- 
scribe the physical conditions in hot DA white dwarf atmospheres correctly. A comparison with previous analyses 
highlights the improvements as well as the limits of our new models. 
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1. Introduction 

White dwarf (WD) atmospheres exhibit a quasi-mono- 
elemental chemical composition. They are either practi- 
cally pure hydrogen or pure helium, with the hydrogen rich 
sequence being called DA and the helium rich sequence be- 
ing summarized as non-DA. In the following, we will re- 
strict our discussion to DA white dwarfs, even though sim- 
ilar physics applies to non-DAs. The basic mechanism for 
pu rification of white dwarf envelopes was first explained 
by Schatzman (1949 ): The strong gravitational field (typ- 
ically, log g — 8 in cgs units) on the surface of WDs yields 
a steep pressure gradient. Pressure driven diffusion sepa- 
rates the elements according to their atomic weight, which 
is generally referred to as gravitational settling or sedi- 
ment atii 



TOTir 



considerable radiative acceleration on heavier elements in 
WD atmospheres, provided the luminosity is high enough 
( IVauclair et al. 1979t [Fontaine fc Michaud 1979| ). In the 
considered temperature range the luminosities are indeed 
still high (i.e. L/Lq > 1). This is sufficient to sustain pho- 
tospheric elements other than hydrogen with abundances 
no higher than 10^"' relative to hydrogen. In optical spec- 
tra of hot white dwarfs, these trace pollutants become di- 
rectly visible only when great effort is put in the detection 



of individual lines (Dupuis et al. 2000). Through resonance 
lines in the UV and even more through the sheer number 
of lines in the EUV, however, these spectral ranges are 
strongly affected. The flux deficit caused by the additional 
opacity in the EUV as compared to pure hydrogen atmo- 



spheres, hinted at already by earlier data (Mewe et al 
1975 ;_ Hcarn et al. 197"6|; Lampton et al. 197(. ; Margon 



range starts trom Jeff ~ 4U, UUUK upwards), however, 
this diffusion process is disturbed by the radiation pres 



In |ftot vvus I, an approximate enec Give rempcracure ^t al. 1976| ; [Shipman 1976D , became more and more ev- 



sure tkat elements experience, depending on the over- 
lap of heir wavelength-dependent opacity with the max- 



imum spectral flux, the overlap ot the maximal spcc- 
tral intensity with the region ot maximal opacity yields a 



idcnt as the quality of observations improved wit h the 
HEAP (Einstein), EXOSAT and ROSAT satell i tes (|Kahn 
et al. 1984|; [Petre et al. 1986|; [Jordan et al. 19871; [Paerels fc 
Heise 1989[ ; [Barstow et al. 1993[ ; [lordan et al. 1994 IWolfl 



et al. 1996[ ). While one EXOSAT ( [Vennes et al. 1989[) and 
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various ROSAT observations revealed the opacity to be 
mainly due to absorbers other than helium, observations 
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with EUVE made possible a more detailed investigation 
of the nature of absorbers. 

The most recent analyses of hot white dwarf EUVE 



spectra have been presented by Barstow et al. (1997) and 
Wolff et al. (199^ ). Using the latest non-LTE and LTE at- 
mosphere models, both groups derived a metal mix scal- 
ing factor for each object in their respective WD sam- 
ple. Barstow et al. use this single scaling factor per object 
to adjust metal abundances predicted for its com bination 
of effe ctive te mperature and surface gravity by [Chayer 



et al. (19954 0), whereas ^Volff et al. use relative metal 



abundances which they have derived for the standard star 
G 191-B2B as a typical metal mix and give the appro- 
priate scaling factor that they call metallicity for each of 
their sample objects. 

In the models used, the distribution of the metals 
was assumed to be homogeneous, i.e. abundances were 
treated as being the same in all atmospheric depths. But 
as the effects of sedimentation and radiative levitation 
lead to a chemical gradient, fitting the emergent flux of 
real, stratified WD atmospheres with homogeneous atmo- 
sphere models has its limits, and true consistency cannot 
be achieved. Therefore, attempts were made to include the 
impact of a chemical gradient via an ad hoc stratification: 
For iron in G 191-B2B ( [Barstow et al. 1999D , which could 
successfully reproduce the EUVE spectrum of G 191-B2B, 
and for nitrogen in RE J1032+535 ( [Holbcrg et al. 19991 ). 

As discussed above, the basic physical mechanisms 
leading to the stratification have long been known, so 
their actual modelling appears preferable to an ad hoc 
approach. Chayer et al. (1995a, ^ were the first to cal- 
culate an extensive model grid for white dwarfs predict- 
ing the chemical stratification assuming equilibrium be- 
tween gravitational settling and radiative levitation. An 
improvement has now been achieved through the intro- 
duction of self-consistent non-NLTE model atmospheres 
which account for the coupling between chemical stratifi- 
cation and the radiation field ( Drcizlcr 1999 ). The success- 
ful applicati on of these models to t he EUVE spectrum of 
G 191-B2B dDrcizler fc Wolff 1999[) has motivated a more 



systematic analysis of a larger sample of hot DA white 
dwarfs which we will present here. 

The first paper on this investigation exploits EUVE 
spectra, further papers will deal with the UV and opti- 
cal wavelength ranges and give detailed metal abundance 
pattern predictions. The availability of EUVE spectra de- 
fines, as in Wolff et al. (1998| ), the overall sample, which is 
more thoroughly introduced in Sect. ^j. The observations 
and particular considerations concerning the interstellar 
medium are shortly described there as well. Sect. || out- 
lines the principle of the model calculations and focuses on 
the characteristics of the model grid computed for com- 
parison with the sample on hand. The methods and re- 
sults from fitting the new class of theoretical spectra to 
the observations are presented in Sect. |[ What implica- 
tions these results involve, and why further efforts in this 
direction are both justified by the present success and at 



the same time necessary due to remaining puzzles, is being 
discussed in Sect. H. 



2. EUVE: Observations of hot DA white dwarfs 

Generally speaking, the maximum spectral flux of hot 
WDs lies in the extreme ultraviolet (EUV) . The fact that 
the main opacity of heavy elements at these temperatures 
tends to lie in the same range means that radiative acceler- 
ation can sustain traces of them in the atmosphere, but it 
also implies that the EUV is the best suited spectral range 
for the detection of these elements. Consequently, we start 
by surveying EUVE spectra, even though its instrument's 
spectral resolutions are not high enough to identify indi- 
vidual lines. Identification element by element is, however, 
not a necessary requirement with the new models: They 
impose no need for the use of a pre-defined metal mix, 
nor do they require the adjustment of e.g. a scaling fac- 
tor, as the abundances are not free parameters but are 
derived from an equilibrium condition (see Sect. As 
the EUV flux is very sensitive to the chemical composi- 
tion of the photosphere, a comparison in this regime is 
extremely useful to test the predictions of the diffusion 
models. Interpreting the repeated analysis as a test case 
for the models, it is reasonable to start out by following 
previous investigations as closely as possible. For this rea- 



son, w e chose to use the 26 DAs already analyzed by Wolfl 
19991) . 



2.1. The EUV selected DA sample 

This sample comprises all hot {Tcs > 40, 000 K) DA white 
dwarfs with EUVE observations, which are mostly taken 
from the EUVE public archive (now available at the mul- 
timission data archive at STScI), and a few from own ob- 
servations. The data sets and the reduction procedures 
applied are described elsewhere ( Wolff et al. 1998 ; Wolfl 
1999); those publications also show the reduced spectra in 
the form used here, i.e. flux calibrated SW, MW and LW 
data combined to yield a single spectrum. 

The effective temperatures of the sample stars lie 
above Tcs — 40, 000 K, where photospheric metal abun- 
dances still are large enough to be detectable at all. Above 
Tcff = 70, 000 K, mass loss effects might start to disturb 
the expected equilibrium between gravitational and ra- 



diative acceleration (Unglaub & Bues 1998), but the sam- 



ple contains no stars hotter than that. Wolff et al. (1998) 



have grouped the objects into four different categories. 
The first set consists of G 191-B2B like objects, where 
the high metallicity yields a steep flux drop towards short 
wavelengths, beyond the interstellar He ii absorption edge. 
The second group comprehends GD 246 and similar ob- 
jects which contain less photospheric metals, but still more 
than those in the third group whose spectra Wolff et al. 
could fit with pure hydrogen atmospheres. A fourth rather 
inhomogeneous set compiles the remaining objects. For de- 
tails, especially on the effective temperature, surface grav- 
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ity and metallicity region in parameter space covered by 
the different groups, we refer once more to [Wolff (1999| ). 



2.2. Treatment of the ISM 

As briefly remarked on before, the interstellar medium can 
considerably attenuate the stellar flux. For the EUV, the 
most important features are the Hi, Hei and Hen bound- 
free ground state absorptions with edges at 9 11. 7 A, 
504.3 A and 227.8 A, respectively. Again, this is bein g ac- 
counted for in the same way as in Wolff et al. (1998 ): The 
effect of given interstellar column densities o n the theoret- 
ical flu xes is modelled and applied following Rumph et al 
(1994 ). The derivation of Hi column density depends on 



the effective temperature of the model used, and the He i 
and Hell column densities on its absorber content. Due to 
the interweavement of Tcfi and N(Hi) and the fact that 
the Lyman edge li es outside the range observed by EUVE, 
Wolff et al. (1998 ) have started from optically determined 



effective temperatures (by Finley et al. 1997), and addi- 
tionally have retained surface gravities from that analysis 
unchanged. In the present approach, we start from the ef- 
fective temperatures from Wolff (1999) and adopt his Hi 
column densities unalteredly whenever possible. 

Even though the He i and He ii column densities have 
to be constantly modified during the fitting procedure de- 
pending on the assumed photospheric parameters, all in- 
terstellar contributions will in the following be treated as a 
secondary effect on the stellar light that can be dealt with 
independently of the analysis of the source itself, disre- 
garding the fact that an appropriate correction can only 
be made after a set of photospheric parameters has been 
adopted. As the interstellar absorption is more impor- 
tant at longer wavelengths while the absorption by photo- 
spheric metals increases towards shorter wavelengths, the 
different contributions can partly be separated. 



3. Equilibrium ansatz: Chemically stratified model 
atmospheres 

Under the assumptions that no processes competing with 
diffusion and levitation are present and that diffusion 
timescales are short, the distribution of trace elements in 
a stellar atmosphere can be expected to take on an equi- 
librium state, as proposed by Chayer et al. (19954 |b|)- In 
such a situation, mean diffusion velocities would be zero 
at any location in the atmosphere, due to an exact bal- 
ance of gravitational, radiative and electrical forces (when 
neglecting thermal and concentration diffusion). A more 
detailed formulation of this condition for polluting ele- 
ments (i) in a hydrogen (or helium) plasma (1) can be 
found elsewhere (Dreizler 199!;). It yields, in a short out- 
line, 



{Ai - Ai)mpg - (Z, 

— AiTOp(7iad,i 

0, 



Zi)eE 



(1) 



where rriiFi are the forces acting on element (i), Ai are 
the atomic weights, Zi are the mean electrical charges, E 
is the electrical field caused by the separation of electrons 
and ions, nip is the proton mass, and f/iad.i is the radiative 
acceleration acting on element (i). Using hydrostatic equi- 
librium and charge conservation, this leads to the simple 
term 



(2) 



One immediately recognizes that requiring diffusion veloc- 
ities to be zero is indeed the same as asking for the effects 
of gravitational settling and radiative levitation to cancel 
each other. Knowing that the above quantities stand for 



^ A,{Zi+i) I y 



(3) 
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where pi is the mass density of the element (i), the 
frequency dependent absorption coefficient which includes 
all contributions of this element at the frequency and 

is the Eddington flux, it furthermore becomes clear 
that (7rad,i = ffcff,i Can only be fulfilled by fixed values for 
Pi. This equilibrium solution for pi is depth dependent. 
For a plane parallel model atmosphere in hydrostatic equi- 
librium this means that equilibrium abundances are solely 
determined by the photospheric parameters Tog and logg, 
in other words, the photospheric abundances are no longer 
free parameters. Calculating the equilibrium abundances 
inevitably leads to a chemically stratified atmosphere. 

Chayer et al. (1995a,^) have presented extensive re- 
sults of such calculations for white dwarfs. Unfortunately, 
their predicted equilibrium abundances were not quantita- 
tively consistent with observations. This might have been 
partly due to the fact that the observations were analyzed 
with chemically homogeneous model atmospheres, mean- 
ing that homogeneous abundances were compared to pre- 
dictions for stratified atmospheres. It is more probable yet 
an effect that arises because the feedback of the modified 
abundances on the radiation field was not being accounted 
for. 

A self-consistent solution of the equilibrium condition 
and the atmosp heric structure has recently been presented 
( Dreizler 1999| ) . In addition to this newly introduced cou- 
pling, it is also a novelty that the calculations are be- 
ing performed under non-L TE conditions. This has turned 
out t o be a crucial point ( Schuh 200C ; Dreizler fc Schuh 
^OOll) . The code itself rehes on an iteration scheme that 
alternates between the determination of new equilibrium 
abundances and the corresponding solution for the atmo- 
spheric structure, which is described in Werner & Dreizlei 
19991) . 



3.1. Scope of application 

It remains to be justified why the white dwarfs in our 
sample will presumably obey the presented equilibrium 
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condition. Processes eligible to disturb or prevent its ad- 
justment can be mass loss, accretion from the interstel- 
lar medium, convection or mixing through rotation. As 
explained by Unglaub & Bues (1998), mass loss would 



have the effect of homogenizing a chemical stratification, 
but as has also been shown in the same paper, mass loss 
rates drop below the critical limit (10~^^ Mo/yr) for DAs 
cooler than Toff ~ 70, 000 K so that this phenomenon 
sh ould not occur in a ny of the sample stars. Calculations 
of MacDonald (1992 ) that treat the interaction between 
accretion and the white dwarf wind reveal that WDs in 
our sample are not affected by accretion. It is prevented 
since Lwd/^o > 1, so that at least a critical mass loss 
rate of 3 • 10"^* to 10"^^ Mo/yr is sustained. 

Convection would of course homogenize abundances, 
too. However, convective instability in the outer layers oc- 
curs only below T^s = 12, COOK for DAs, i.e. well below 
Teff — 40, 000 K. Rotation could lead to a mixing through 
meridional currents; but as most white dwarfs are very 
slow rotators (Heber et al. 1997; Koester et al. 1998), this 
is not likely to interfere cither. 

Even if it is left undisturbed, it takes an atmosphere 
a certain time to reach an equilibrium state. Due to the 
high surface gravities, diffusion time scales are of the or 
der of months in the outer layers of white dwarfs (Koester 
1989|P)ased on diffusion coefficients by Paquette et al 



1986^ This is nothing but an instant compared to evo- 
lution time scales (~ 10^ yr), and this is the reason why 
DA atmospheres in the effective temperature range be- 
tween 40, 000 K and 70, 000 K can be regarded to have 
their photospheric abundances set to equilibrium values 
at any time of their evolutionary stage. In the temper- 
ature range above 40, 000 K the flux maximum still lies 
in the EUV, i.e. in the same range where multiple lines 
of heavier elements provide substantial opacity, making 
radiative levitation an efficient mechanism for sustaining 
absorbers. It should in this context also be noted that the 
evaluation of the equilibrium abundances does not include 
the contributions of thermal and concentration diffusion, 
which have much longer time scales. 



3.2. The model grid 

The model grid spans the Tcff and logg plane as sug- 
gested by the results from 



Wolff (1999) 



The effective 

temperatures cover a range from 38, 000 — 70, 000 K in 
2,000 — 3, OOOK steps, the logarithmic surface gravity a 
range from 7.2 — 8.3 in 0.1 — 0.2 dex steps. The grid is 
however not complete and concentrates around the pre- 
viously found parameter combinations for the objects as 
listed in Table |^. We emphasize once more that the total 
absorber content, i.e. the entire depth-dependent distri- 
bution of each individual trace element, is solely defined 
by a models' [TefT,logf/] - combination, which makes Tcfi 
and log 5 the only free parameters. The elements included 
in the model calcu lations are the same as those used by 
Wolff et al. (1998 ) for the definition of their metal mix. 



which in turn are those that they were able to identify in 
HST GHRS data of the standard G 191-B2B. Other ele- 
ments were assumed to be present in abundances too low 
to significantly contribute to the EUV opacity (however, 
stratification may well invalidate this assumption). While 
other elements than those used in Wolff's and this analy- 
sis have been identified in UV spectra of some of the sam- 
ple stars, including G 191-B2B itself, the low abundances 
of 2.5 • 10^^ for phosphorus and 3.2 • 10^^ for sulfur, for 
examp le, as reported from the analysis of ORFEUS spec- 
tra by Venues et al. (1996a ), confirm that these element's 
opacities will presumably not contribute significantly to 
the total EUV opacity. To investigate the possible effects 
of several additional elements anyhow, (which currently 
results in a dramatic increase of CPU time demand) , fur- 
ther analyses will have to be conducted (see also Sect. 5.2). 
For this analysis, we aimed at being able to compare our 
results to a previous analysis in the EUV, as mentioned 
before. Consequently, the list of elements considered be- 
sides H reads He, C, N, O, Si, Fe and Ni. Differences 
remain with regard to the number of atomic levels and 
with regard to LTE/non-LTE population of these levels. 
Details for the model calculations presented here are sum- 
marized in Table|^. Note that so far, the model atoms used 
for the solution of the atmospheric structure and for the 
solution of the equilibrium condition are identical, which 
implies that the relatively high accuracy required for the 
correct evaluation of the radiative acceleration slows down 
the atmospheric structure calculation. 

The time it takes a model with the specified atomic 
data to converge depends primarily on the stratification 
of the input model. For homogeneous input models with 
abundances derived from the published mctallicities as 
a first guess, roughly 5x10^ CRAY-CPU seconds are re- 
quired, while using converged stratified models with pho- 
tospheric parameters that are somewhat off the desired 
ones yields convergence after of the order of 5x10^ CPU 
seconds. Due to this particularly large CPU time demand 
per model, the convergence criterion for the models was 
set to be a less than 5% change in the model flux, which 
means that the maximal relative correction at any fre- 
quency point is no larger than that limit, compared to 
the result from the previous diffusion iteration step. Only 
single lines (usually two or three He lines) are actually af- 
fected by those larger variations, the continuum flux itself 
is much more exactly determined and easily complies to 
our usual "< 10"''" criterion. 

The next generation of models will benefit from ma- 
jor improvements of the model atmosphere code and the 
iteration scheme for the equilibrium condition. Both will 
reduce the CPU time drastically. 

As to the uniqueness of a particular solution, tests in- 
dicate that the results agree to within our error limits 
regardless of the extremely different stratification of di- 
verse start models, making a strong case at least for the 
stability of the result. We therefore trust that, no matter 
what the initial conditions in the start model are, the fi- 
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Table 1. Summary of the model atoms. Super-levels 
(marked by *) comprise hundreds to thousands of atomic 
levels. Line transitions between these super-levels (marked 
by "*") are composed of all transitions between the indi- 
vidual levels comprising the super-levels (see Werner & 
Dreizler 1999 for details). 
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He 
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O 


Si 


Fe 


Ni 


non-LTE 


















levels 


17 


54 


70 


50 


24 


25 


28* 


29* 


lines 


56 


106 


218 


114 


35 


53 


77+ 


44+ 



tation of depth dependent stratifications for all elements 
included in the calculations will be published separately. 

4. Model matching: Testing the new models and 
re-analyzing the DA sample 

To directly compare a theoretical flux distribution with 
the observed one, it is first re-binned in wavelength and 
normalized to the observed visual magnitude. T he values 
for my have b een chosen using the compilation of McCook 
|fc Sion (1999|) and are listed with their individual refer- 
ences in Table |[ Then the effect of in terstellar absorption 
is calculated as sketched in Sect. 2.2. 



nal abundances in the converged model will always be the 
same. 

The emergent fluxes used for comparison with the ob- 
servations have been calculated from the model atmo- 
spheres on a frequency grid optimized for the EUV spec- 
tral range. This grid covers A = 50 - 700 A in detail and 
contains more than 30,000 frequency points in total. 

3.3. Characteristics of the models 

The models show the general properties as expected 
from diffusion theory, i.e. the overall trace element abun- 
dances decrease with lower T^s and higher \ogg values. 
Accordingly, the emergent spectra approach the flux dis- 
tributions of those of pure hydrogen model atmospheres 
at sufficiently evolved parameters on the cooling sequence. 
The vertical element distribution however does not gener- 
ally follow the simple picture of monotonously increasing 
local abundances with depth, but is strongly affected by 



the respective radiative acceleration (see Dreizler 1999, 
Schuh 2000, or Dreizler & Schuh 2001 for exemplary 
graphical representations). In particular, gradients may 
change signs repeatedly, meaning this can be brought 
about by diffusion processes alone, in contrast to other 
notions that additionally evoke mass loss effects to ex- 
plain local abundances that increase towards outer at- 
mospheric regions (e.g. [Holberg et al. 1999| ). In a nut- 
shell, even equilibrium calculations yield rather non-trivial 
abundance distributions. 



As in p revious calculations ( Dreizler 1999 ; Dreizler fc 
Wolff 1999| ), the iron group element abundance can reach 
solar values or more, with nickel being at virtually the 
same level as the iron abundances. Though their domi- 
nant effect on the EUV opacity is partly due to the many 
lines these elements exhibit, it is equally important that 
their abundances are, over a wide parameter range, simply 
larger than those of other metals. Of course these two ef- 
fects are tightly correlated, since more lines directly trans- 
late into a stronger radiative acceleration, which in this 
case easily compensates for the higher atomic weights and 
thus a stronger response to gravity. 

This represents just a short overview of some general 
effects that can be seen in the models; detailed presen- 



4.1. Comparison of theoretical and observed spectra 

The model matching was performed on the mentioned grid 
which does not in all parameter ranges correspond to a 
tight mesh, and no interpolation was used. The quality of 
an individual fit was evaluated by eye. 

Of course, changing the interstellar He column den- 
sities seemingly has a big impact on the overall run of 
the spectrum. However, the region most sensitive to metal 
abundances lies at least partly beyond the Heii absorp- 
tion edge, so once this and the lower energy edges are 
fitted all remaining deviations have to be due to either a 
different flux level as caused by the effective temperature 
or a different amount of absorber content as regulated by 
the surface gravity in the equilibrium models. 



4.2. Revised atmospheric parameters 

The present model matching led to new atmospheric pa- 
rameters where mainly the surface gravities differ from 
former specifications. The improvement of many of the 
individual fits (see Fig. p|) as compar e d to fits wi t h homo- 
geneous mode ls (in Wolff et al. 199S ; Wolff 199£ ; Dreizler 
|fc Wolff 199£ ) vindicates these adjustments. 

With homogeneous models, discrepancies were encoun- 
tered especially for higher metallicity stars in the shorter 
wavelength regimes, which could not be adjusted by a dif- 
ferent metallicity value as this would have deteriorated 
the quality of the fit at other wavelengths. This difficulty 
is resolved with the new models, as they offer different 
abundances at different formation depths. Consequently 
they can often better reproduce especially these formerly 
delicate regions. Fig. ^ shows an example. 

Fig. |2| displays all observations with the best-fit model 
spectra over-plotted, both shown at lA resolution. Table|| 
lists the parameters of the best-fit model for each object, 
along with the interstellar column densities used to pro- 
duce the plots. These are to be regarded as preliminary 
as fitting improved models to the spectra may yet again 
yield different values. High Hen column densities, found 
throughout the sample, may be an indication that some 
opacities are still missing. It additionally lists a quantity 
denoted mi, which stands for the metal index that is to 
be introduced now. 
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Fig. 1. EUVE spectrum of MCT 2331-4731 with the best-fit model at Tcff = 56, 000 K and logg = 7.7 overlaid 
(dashed). The photon flux is displayed on a logarithmic scale. 
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Table 2. New and previous results 
is given in lO^^cm"^. See Section 4, 



( Wolff et al. 199^ ) of the EUV analysis ordered by decreasing metal index mi. Nh 
for an explanation of the theoretical metal index mi. 



WD-N° 


Name 


mv 




log g 


Tcff 


mi 
new 


Nh 
results 


Hel/H 


Hell/H 


log g Tcff m 
previous results 


WD 0621-376 


REJ0623-377 


12.089 


i 


7.2 


61000 


3.49 


5.0 


0.09 


0.15 


7.27 


58000 


2.0 


WD 2211-495 


REJ2214-493 


11.7 


1 


7.4 


66000 


3.02 


5.8 


0.07 


0.15 


7.38 


66000 


4.0 


WD 0232+035 


Feige 24 


12.56 


2 


7.6 


59000 


1.22 


2.72 


0.068 


0.25 


7.17 


58000 


1.0 


WD 0455-282 


MCT 0455-2812 


13.95 


3 


7.8 


66000 


1.20 


1.3 


0.063 


0.3 


7.77 


66000 


1.0 


WD 0501+527 


G 191-B2B 


11.79 


2 


7.6 


56000 


0.99 


2.05 


0.071 


0.3 


7.59 


56000 


1.0 


WD 2331-475 


MCT 2331-4731 


13.1 


4 


7.6 


56000 


0.99 


8.5 


0.08 


0.1 


8.07 


56000 


0.75-1. 


WD 1056+516 


LB 1919 


16.8 


4 


8.2 


70000 


0.61 


16.0 


0.04 


0.05 




69000 


0.1 


WD 1123+189 


PC 1123+189 


14.13 


1 


7.9 


54000 


0.43 


11.9 


0.09 


0.052 


7.63 


54000 


0.4 


WD 0027-636 


MCT 0027-6341 


15.0 


4 


8.2 


64000 


0.42 


21.5 


0.068 


0.06 


7.96 


64000 


0.2 


WD 1234+482 


PC 1234+482 


14.38 


5 


8.1 


56000 


0.31 


11.7 


0.09 


0.05 


7.67 


56000 


0.2 


WD 2309+105 


CD 246 


13.09 


2 


8.2 


56000 


0.25 


18.0 


0.05 


0.03 


7.81 


59000 


0.25 


WD 0131-164 


CD 984 


13.98 


6 


8.1 


50000 


0.20 


22.0 


0.068 


0.035 


7.67 


50000 


0.2 


WD 2247+583 


Lanning 23 


14.26 


7 


8.3 


56000 


0.20 


40.0 


0.068 


0.03 


7.84 


59000 


0.25 


WD 0004+330 


CD 2 


13.85 


2 


8.25 


50000 


0.14 


82.4 


0.068 


0.01 


7.63 


49000 


< 0.1 


WD 2111+498 


CD 394 


13.09 


8 


7.9 


40000 


0.13 


6.5 


0.07 


0.15 


7.94 


39600 


0.25 


WD 1314+293 


HZ43A 


12.914 


9 


8.3 


50000 


0.13 


0.9 


0.06 


0.0 


7.99 


50800 


0.0 


WD 1543-366 


RE J 1546 -364 


15.81 


10 


8.3 


50000 


0.13 


50.4 


0.068 


0.03 


8.88 


45200 


0.0 


WD 1057+719 


PC 1057+719 


14.68 


1 


8.0 


42000 


0.12 


20.7 


0.068 


0.052 


7.90 


41500 


0.1 


WD 0630-050 


REJ0632-050 


15.537 


1 


8.2 


44000 


0.09 


30.1 


0.01 


0.055 


8.39 


44100 


0.0 


WD 1631+781 


REJ1629+780 


13.03 


11 


8.15 


42000 


0.09 


35.0 


0.068 


0.0 


7.79 


44600 


0.05 


WD 0838+035 


REJ0841+032 


14.48 


1 


8.1 


40000 


0.08 


16.0 


0.068 


0.052 


7.78 


38400 


0.0 


WD 2321-549 


REJ2324-544 


15.2 


1 


8.2 


42000 


0.08 


9.5 


0.05 


0.06 


7.94 


45000 


0.05 


WD 1029+537 


REJ1032+535 


14.455 


1 


8.3 


44000 


0.08 


7.5 


0.005 


0.04 


7.77 


44000 


0.0 


WD 2152-548 


REJ2156-546 


14.44 


1 


8.3 


44000 


0.08 


7.0 


0.04 


0.04 


7.91 


44000 


0.0 


WD 0715-703 


REJ0715-705 


14.178 


1 


8.3 


44000 


0.08 


21.9 


0.068 


0.015 


8.05 


44000 


0.0 


WD 2004-605 


REJ2009-605 


13.6 


1 


8.3 


42000 


0.06 


17.5 


0.06 


0.015 


8.16 


41900 


0.0 



^ iMarsh et al. (1997) ^ 


Kidder et al. (1991) ^ 


Barstow et al. (1994) Pounds et al. (1993) ^ 


Green (198C) ^ Wesemael (1995) 


^ [Vennes et al. (1997 




Bergeron et al. (1992) ^ Bohlin et al. (1995 


) Vennes et al. (1996t) Schwartz et al. (1995) 
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Fig. 2. EUVE spectra of the program stars with theoretical spectra overlaid (dashed), ordered by decreasing metal 
index mi. We only show the spectral range with significant signal. Parameters of the models as well as interstellar 
column densities for hydrogen and helium can be found in TableH. 
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4.3. Metal index 



5. Discussion 



With depth dependent abundances, it is difficult to quote 
one representative value for each element. To evade this 
problem, we have introduced a photospheric parameter de- 
pendent quantity that relies entirely on the predictions by 
diffusion theory. As explained above, the absorber content 
of the models is determined from the equilibrium between 
the radiative and the gravitational forces. The radiative 
acceleration scales with T^g due to the occurrence of the 
Eddington flux in the lower expression of Eq. ||, the effec- 
tive gravitational acceleration scales with g as is evident 



Our new models can reproduce the observed EUV spectra 
of hot DA white dwarfs. Since the new models predict the 
chemical composition from the equilibrium between sedi- 
mentation and radiative acceleration, the number of free 
parameters is drastically reduced to the effective temper- 
ature and surface gravity only. The good agreement is a 
strong evidence that the interplay of these two processes 
defines the chemical composition and stratification. The 
agreement is, however, not in all cases better than with 
chemically homogeneous models, which requires a more 
detailed discussion. 



from tlic upper expression in the same equation. The ab- et al. 



Stars sim ilar to G 191-B2B (the group one from Wolfl 



sorber content should thus be comparable along lines of 
constant T^g/g. To map this ratio onto a dimensionless 
and easier-to-read parameter, we define 



mi 



4.10-i2.T,V5 /[KV/cm] 



(4) 



where the leading factor has been chosen to yield a result 
of the order of one for the photospheric parameters of 
the standard star G 191-B2B. Lines of constant mi are 
displayed in Fig.^. 



It should be emphasized that this quantity is not de- 
rived from the actual calculated equilibrium abundances 
in the models. It is based on a much simpler evaluation of 
just the same idea that underlies the construction of the 
models. Knowing this, it is all the more surprising how 
well the metal indices mi at the newly determined pa- 
rameters match the metallicities m that Wolff (199^ ) has 
derived for the objects. 



19981) can be reproduced significantly better with 
our self-consistent stratified models. This result is not sur- 
prising since the first and successful application of these 
models to G 191-B2B (|Dreizler fc Wolff 1999[) motivated 



this work in the first place. Our new models can also re- 
prod uce the stars classi fied as pure hydrogen (group three 
from Wolff et al. 1998 ), demonstrating that the vanish- 



ing of trace amounts of metals is correctly reproduced. 
Within this group, however, the surface gravity had to be 
increased by up to 0.5 dex in order to achieve a good fit 
(see also Table |). 

The fits of PG 1123+189, GD246 and GD984 are 
not completely satisfying but deviations are of the same 
order as compared to fits with homogeneous models. 
The fits for HZ43A, RE J1032-H535, REJ2156-546 and 
REJ2009— 605 are not quite as good as the ones which 
could be achieved with spectra derived from the grid of 
homogeneou s mo dels. Further improvements of our mod- 
els (see Sect. 5^) may resolve these slight discrepancies. 

In four cases (LB 1919, MCT 0027-6341, 
PG 1234+482, and GD394), using the new models 
results in a worse fit than before. One may be tempted 
to explain this by the fact that our model grid is not 
yet extended enough to find a reasonable agreement. 
However, significant individual deviations between models 
and observation could also originate due to physical rea- 
sons. As discussed in Sect. 3.1, competing processes can 



disturb the equilibrium between radiative acceleration 
and gravitational settling. As is evident from the overall 
good fit of our models, these are unimportant in general 
but individual exceptions are possible. In the case of 
LB 1919 for example, an unusually high rotation rate 
has been suspected (Finlcy et al. 1997). Prior to a final 
decision, we will have to implement several improvements 
to our models. 



5.1. Remarks on g and Qrad 

As mentioned above, our surface gravities are in several 
cases higher than in earlier analyses. Those analyses had 
evaluated the Balmer lines for a determination of the sur- 
face gravity. Should the surface gravities determined with 
our models prove to be reliable, they would open up a new 
possibility: The surface gravity governs the gravitational 



S.L. Schuh et al.: Equilibrium abundances in hot DA white dwarfs I. 



9 




Fig. 3. Revised parameters of the program stars in the [Toff, log g] plane. Lines indicate constant values of the metal 
index as defined in Sect.^^, mi = 4., 0.65, 0.15, 0.05, respectively. The shaded area indicates the metallicity of DA 
white dwarfs of the GD 246 group of Wolff et al. (1998| ). These as well as the objects classified as irregular fall into 
this regime and are marked by squares. Stars above that region belong to the G 191-B2B-group (identified by starry 



symbols), stars below to the pure- hydrogen-group (marked by diamond-shaped symbols). See also Sect. 2.1 



settling. Through the equilibrium condition, the amount 
of trace metals in the atmosphere is determined. In ef- 
fect, we therefore evaluate the metal abundances in order 
to derive the surface gravity. The new models with strati- 
fied metal abundances can modify the Balmer lines via the 
changed atmospheric structure. This could explain the de- 
viations in the surface gravity determination. A more de- 
tailed investigation is in progress. Our proposed method 
would be a very sensitive indicator but it is, unfortunately, 
also sensitive to systematic errors in the calculation of the 
radiative acceleration. The shifts in log g necessary to ob- 
tain good fits for several of the objects consequently could 
indicate that the accuracy of the grad calculation needs to 
be further improved. The next generation of models will 
therefore aim at determining radiative accelerations with 
higher precision (see Sect. 5.2). 



One independent result, however, might be worth con- 
sidering in this context: Detailed parallax and gravita- 
tional redshift measurements with HST of the binary sys- 
tem Feige24 ( [Benedict et al. 2000| ; [Vcnncs et al. 2000[ ) 
also yielded a higher surface gravity for the WD com- 
ponent than several optical analyses (log g=7.2-7.5). The 
consistency of our result (log g—7.7) with Vennes et al.'s 
(log g— 7 .5-7 .7) might be interpreted as a hint that our 
current results do not suffer from an inaccurate grad de- 
termination in the models after all. Instead, analyses with 
homogeneous models would have been biased due to sys- 
tematic errors introduced through the neglect of stratifi- 



cation. But even though this may sound plausible, it is 
too early to call such a statement anything more than a 
suspicion. 

5.2. Conclusions and perspective 

In this paper, we have presented our chemically strati- 
fied model atmospheres describing the chemical compo- 
sition through the equilibrium between gravitational set- 
tling and radiative acceleration self-consistently. The ap- 
plication to 26 EUVE spectra of hot DA white dwarfs has 
revealed an overall good agreement, demonstrating the po- 
tential of these new models. 

In upcoming papers, we will check the derived param- 
eters through a re-analysis of the Balmer lines and of the 
metal lines observable in UV spectra. The analysis of opti- 
cal spectra with our new stratified models should allow to 
check for systematic error as compared to previous anal- 
yses with chemically homogeneous models. This will also 
be an important test for the consistency of atmospheric 
parameters derived from different parts of the spectrum, 
which has often posed a problem in the past, and to which 
we hope the new type of models will ultimately be a so- 
lution. Only this test will show on how a firm ground the 
parameters presented here stand. Likewise, the determina- 
tion of the interstellar column densities has to be regarded 
as preliminary, since these cannot be determined indepen- 
dently from the theoretical spectra. 
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The current models are only a first step, in the future 
we will be able to present an improved set of stratified 
model atmospheres. A "next generation" version of our 
stellar atmosphere code itself as well as a more efficient 
coupling between the computation of the model structure 
and the chemical stratification will allow to include much 
more detailed model atoms and at the same time many 
more chemical elements. This influences and is crucial for 
the accuracy of the calculation of the radiative accelera- 
tion. By evaluating the scale of the deviation in compari- 
son to the current results it should then become clear how 
trustworthy the latter are. 
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